#include <iostream>

#include "Vector/vector.03.h"

#include "icc.01.h"

using namespace std;

int main() {
	
	
	ICC icc;
	
	uint N = 10000;
	double r1 = 0;
	double r2 = 1;
	
	//vector< double > r = icc.getr<double>(N, r1, r2);
	
	/*
	for (uint i=0; i<N; i++)
		printf("% .3e\n", r[i]);
	*/
	/*
	vector <int> bin(100, 0);
	double dr = (r2-r1) /bin.size();
	for (uint i=0; i<r.size(); i++) 
		bin[int((r[i]-r1)/dr)]++;
	
	for (uint i=0; i<bin.size(); i++) {
		double r = i*dr;
		double dV = (4.0/3*3.141592654) * (3*r*r*dr + 3*r*dr*dr + dr*dr*dr);
		printf("%3e %2i %5e\n", r+0.5*dr, bin[i], icc.m*bin[i]/dV);
	}
	*/
	 
	
	
	//vector< Vector<double> > R = icc.getR<double>(N, r1, r2);
	vector< Vector<double> > R = icc.getHernquistR<double>(N, 0.27);
	for (uint i=0; i<N; i++)
		printf("% .3e % .3e % .3e\n", R[i].x, R[i].y, R[i].z);
	
	return 0;
}

